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Isolated singularities for a semilinear equation for the 
fractional Laplacian arising in conformal geometry 

DelaTorre, Azahara* 

Gonzalez, Marfa del Mar^ 


Abstract 

We introduce the study of isolated singularities for a semilinear equation involving 
the fractional Laplacian. In conformal geometry, it is equivalent to the study of 
singular metrics with constant fractional curvature. Our main ideas are: first, to set 
the problem into a natural geometric framework, and second, to perform some kind 
of phase portrait study for this non-local ODE. 

1 Introduction and statement of results 

We consider the problem of finding radial solutions for the fractional Yamabe problem 
in R”, n >2, with an isolated singularity at the origin. This means to look for positive, 
radially symmetric solutions of 

(—A)^r(; = in R" \ {0}, (1.1) 

where Cn,^ is any positive constant that, without loss of generality, will be normalized 
as in Proposition 2.7. Unless we state the contrary, 7 G (0, ^). In geometric terms, 
given the Euclidean metric \dx\'^ on R", we are looking for a conformal metric 

gw = \dx\'^, re > 0, (1.2) 

with positive constant fractional curvature Q®™ = Cn,7, that is radially symmetric and 
has a prescribed singularity at the origin. 
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Because of the well known extension theorem for the fractional Laplacian [12, 
13, 14] we can assert that equation (1.1) for the case 7 G (0,1) is equivalent to the 
boundary reaction problem 

{ -div(j/“VVF) = 0 in 

VF = w on K" \ { 0 }, ^^ 3 ^ 

n + 27' 

—dj lim y’^dyW = Cn on R” \ { 0 }. 

Here the constant d-y is defined in (2.4). We note that it is possible to write W = 
w, where Kj is the Poisson kernel (2.3) for this extension problem. 

It is known that wi{r) = r 2 ^ together with Wi = K-y *x wi, is an explicit 
solution for (1.3); this fact will be proved in Proposition 2.7 and as a consequence we 
will obtain the normalization of the constant Cn,-y- Therefore, wi is the model solution 
for isolated singularities, and it corresponds to the cylindrical metric. 

In the recent paper [10] Caffarelli, Jin, Sire and Xiong characterize all the nonneg¬ 
ative solutions to (1.3). Indeed, let W be any nonnegative solution of (1.3) in 
and suppose that the origin is not a removable singularity. Then, writing r = \x\ for 
the radial variable in M", we must have that 

W{x,t) = W{r,t) and drW{r,t) <0 V0<r<oo. 

In addition, they also provide their asymptotic behavior. More precisely, if w = IT(-, 0) 
denotes the trace of W, then near the origin one must have that 

n —27 n —27 

Cir“ 2 < w{x) < C 2 r~ 2 , (1.4) 

where ci, C 2 are positive constants. 

We remark that if the singularity at the origin is removable, all the entire solutions 
for (1.3) have been completely classified by Jin, Li and Xiong [43] and Chen, Li and 
Ou [18], for instance. In particular, they must be the standard “bubbles” 

n-2~/ 

w{x) =c( ^ -r) , c, A > 0, Xo G K”. (1.5) 

-I- ja; - xoj/ 

In this paper we initiate the study of positive radial solutions for (I.l). It is clear 
from the above that we should look for solutions of the form 

w^r) = r 2 -^ v{r) on R" \ { 0 }, ( 1 . 6 ) 

for some function 0 < ci < u < C 2 . In the classical case 7 = 1, equation (1.1) reduces 
to a standard second order ODE. However, in the fractional case (1.1) becomes a 
fractional order ODE, so classical methods cannot be directly applied here. 

The objective of this paper is two-fold: first, to use the natural interpretation of 
problem ( 1 . 1 ) in conformal geometry in order to obtain information about isolated 
singularities for the operator (—A)''' from the scattering theory definition. And sec¬ 
ond, to take a dynamical system approach to explore how much of the standard ODE 
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study can be generalized to the PDE (1.3). In the forthcoming paper [23], that is a 
joint paper together also with M. del Pino and J. Wei, we conclude this study by using 
a variational method directly to construct solutions to ( 1 . 1 ) that generalize the well 
known Delaunay (sometimes called Fowler) solutions for the scalar curvature problem. 
Both papers complement each other. 

Before stating our results we need to introduce the geometric setting. On a gen¬ 
eral Riemannian manifold (M”,g), the fractional curvature Q® is defined from the 
conformal fractional Laplacian P® as Q® = P®(1), and it is a nonlocal version of the 
scalar curvature (which corresponds to the local case 7 = 1). The conformal fractional 
Laplacian is constructed from scattering theory on the conformal inhnity M" of a con¬ 
formally compact Einstein manifold (X"+^, g+) as a generalized Dirichlet-to-Neumann 
operator for the eigenvalue problem 

-Ag+U-s(n-s)[/ = OinX, s=^+j, (1.7) 

and it is a (non-local) pseudo-differential operator of order 2"f. This construction 
is a natural one from the point of view of the AdS/CFT correspondence in Physics 
([1,65]). 

The main property of P® is its conformal invariance; indeed, for a conformal change 

4 

of metric gw = we have that 

n + 27 

P^{fw), for all / smooth, ( 1 . 8 ) 

which, in particular when / = 1 , reduces to the fractional curvature equation 

n + 27 

P®(w) = . (1.9) 

The fractional Yamabe problem for equation (1.9) on compact manifolds was consid¬ 
ered in [35, 36], while the fractional Nirenberg problem was introduced in [43, 44, 27]. 
In addition, the study of the singular version for the fractional Yamabe problem was 
initiated in [34], where the authors construct model singular solutions. Here we look 
at the case of an isolated singularity. 

We note that, in the Euclidean case, coincides with the standard fractional 

Laplacian (—A)^, and thus, imposing the constant curvature condition in equation 
(1.9) yields our original problem (1.1). Moreover, one can check that the extension 
problem (1.3) comes from the scattering problem (1.7) when we take to be the 
hyperbolic metric g~^ = £llLAhl£L_ gee Section 2 for a review of the known results on 
the fractional conformal Laplacian. 

We present now the natural coordinates for studying isolated singularities of (1.1). 
Let M = ]R"'\{0} and use the Emden-Fowler change of variable r = e“*, t G K; with 
some abuse of the notation, we write v = v(t). Then, in radial coordinates, M may 
be identified with the manifold R x for which the Euclidean metric is written as 

\dx\^ = dr'^ -br-^5s"-i = +5S"-i] =: e“^‘5o- (1-10) 
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Since the metrics \dx\^ and go are conformally related, we prefer to use go as a back¬ 
ground metric and thus any conformal change ( 1 . 2 ) may be rewritten as 

4 2 4 

gv = \dx\ =Z;"- 2 rgo, 


where we have used relation (1.6). Looking at the conformal transformation property 
for P® given in ( 1 . 8 ) and relation ( 1 . 6 ) again, it is clear that 

pso(^) = = r^(-A)^w, (1.11) 

and thus the original problem ( 1 . 1 ) is equivalent to the following one: fixed go as a 
background metric on R x S" find a conformal metric = u^-^-rgo of positive 
constant fractional curvature Q®”, i.e., find a positive smooth solution v for 

p3o(i;) = on RxS"-b ( 1 . 12 ) 

The point of view of this paper is to consider problem (1.12) instead of (1.1), since 
it allows for a simpler analysis. In our first theorem we compute the principal symbol 
of the operator P^o qj^ x using the spherical harmonic decomposition for S"“^. 
With some abuse of notation, let = —k{k -f n — 2), fc = 0,1, 2,... be the eigenvalues 
of Agn-i, repeated according to multiplicity. Then, any function on R x may be 
decomposed as J2k''^k{t)Ek, where {Ek} is a basis of eigenfunctions. We show that 
the operator P^° diagonalizes under such eigenspace decomposition, and moreover, it 
is possible to calculate the Fourier symbol of each projection. Let 

w(0 =/ e-^^'*v{t)dt (1.13) 

V 27r Jr 

be our normalization for the one-dimensional Fourier transform. 


Theorem 1.1. Fix 7 € (0, ^) and let P^ be the projection of the operator P®“ over 
each eigenspace {Ek). Then 

lf^k) = Q^,{(,)vk, 

and this Fourier symbol is given by 


0^(0 = 2 "^ 


rU + i 


'J ^2 I { 
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F i - 2 
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(1.14) 


Since we are mainly interested in radial solutions v = v{t), in many computations 
we will just need to consider the symbol for the first eigenspace k = 0 (that corresponds 
to the constant eigenfunctions): 
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Now we look at the question of finding smooth solutions u = u(t), 0 < Ci < 
V < C2, for equation (1.12), and we expect to have periodic solutions. The local case 
7 = 1 , presented below, provides some motivation for this statement and as we have 
mentioned, in the forthcoming paper [23] we construct such periodic solutions from 
the variational point of view. Here we look at the geometrical interpretation of such 
solutions and provide a dynamical system approach for the problem. 

First, note that 

—Aw = Cn,i'w^^, w > 0, (1.15) 

is the constant scalar curvature equation for the metric g^- It is well known ([9, 47]) 
that positive solutions of equation (1.15) in ]R"\{0} must be radially symmetric and, 
if the singularity at the origin is not removable, then the solution must behave as in 
(1.4). After the change (1.6), our equation (1.15) reduces to a standard second order 
ODE for the function v = v{t): 


■■ I (n-2)^ ("-2)^ ^ n /"I 1C^ 

— u + ^ V = ^ , u > 0. (1.16) 

This equation it is easily integrated and the analysis of its phase portrait gives that 
all bounded solutions must be periodic (see, for instance, the lecture notes [61]). More 
precisely, the Hamiltonian 

+ (L17) 

is preserved along trajectories. Thus looking at its level sets we conclude that there 
exists a family of periodic solutions {vl} of periods L € (Li, oo). Here 



is the minimal period and it can be calculated from the linearization at the equilibrium 
solution Vl = 1. These {r’i} are known as the Fowler ([30]) or Delaunay solutions for 
the scalar curvature. 

Delaunay solutions are, originally, rotationally symmetric surfaces with constant 
mean curvature and they have been known for a long time ([24, 25]). In addition, let 
E € be a noncompact embedded constant mean curvature surfaces with k ends. It 
is known that any of such ends must be asymptotic to one of the Delaunay surfaces 
([46, 55]), which is very similar to what happens in the constant scalar curvature 
setting (see, for instance [47]), where any positive solution of the constant scalar 
curvature equation (1.15) must be asymptotic to a precise deformation of one of the 
Vl- Delaunay-type solutions have been also shown to exist for the first eigenvalue in 
Serrin’s problem [62]. 

Let us comment here that the constant mean curvature case is very related to 
the constant Q 1/2 case. This is because, fixed a compact manifold with 

boundary M" such that there exists a defining function p so that 5 + = g /is asymp¬ 
totically hyperbolic and has constant scalar curvature Rg+ = —n{n + 1), then the 
conformal fractional Laplacian on M with respect to the metric g = g\M is given by 
^ 1/2 ~ ^ where ^ is the Neumann derivative for the harmonic extension 


5 







and T-Lg the mean curvature of the boundary ([41, 14]). Thus coincides in this 
particular setting with T-Lg up to a multiplicative constant. However, there is a further 
restriction since in the present paper we consider only rotationally symmetric metrics 
on the boundary and thus, not allowing full generality for the original constant mean 
curvature problem. 


Going back to (1.16), we would like to understand how much of this picture is 
preserved in the non-local case, so we look for radial solutions of equation (1.1), which 
becomes a fractional order ODE. On the one hand, we formulate the problem through 
the extension (1.3). This point of view has the advantage that the new equation is 
local (and degenerate elliptic) but, on the other hand, it is a PDE with a non-linear 
boundary condition. Note that because we will be using the extension from Theorem 
2.4 for the calculation of the fractional Laplacian, we need to restrict ourselves to 
7 € (0,1) at this stage. 

The first difficulty we encounter with our approach is how to write the origi¬ 
nal extension equation (1.3) in a natural way after the change of variables r = e~*. 
Looking at the construction of the conformal fractional Laplacian from the scattering 
equation (1.7) on a conformally compact Einstein manifold (X"+^,(;+), we need to 
find a parametrization of hyperbolic space in such a way that its conformal infinity 
M" = {p = 0} is precisely (ffi. x S"“^,po)- The precise metric on the extension is 

9^ = glP^ for 

g = dp^ (1-19) 

where p G (0, 2) and t G M. The motivation for this change of metric will be made 
clear in Section 3. In the language of Physics, g'^ is the Riemannian version of AdS 
space, a model solution of Einstein’s equations which is important in the setting of 
the AdS/CET correspondence since AdS space-time is the right background when 
studying thermodynamically stable black holes [45, 65]. 

Rewriting the equations in this new metric, our original equation (1.3), written in 
terms of the change (1.6), is equivalent to the extension problem 


- divgip^^VgV) + E(p)V = 0 in (A"+i, p), 

V = V on {p = 0}, 

n + 2~f 

—dj lim p'^dpV = on {p = 0}, 

p —>-0 


( 1 . 20 ) 


where the expression for the lower order term E{p) will be given in (5.2). We look for 
solutions V to (1.20) that only depend on p and t, and that are bounded between two 
positive constants. 

We show first that equation (1.20) exhibits a Hamiltonian quantity that generalizes 
(1.17). 

Theorem 1.2. Fix 7 G (0,1). Let V be a solution of (1.20) only depending on t and 
p. Then the Hamiltonian quantity 

H^{t) - ^ {ei{p){dtVf - e{p){dpVf - e 2 {p)V^] + ( 1 . 21 ) 
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is constant with respect to t. Here we write 


eiip) = (l - (n - 2 + nli) , 

and the constant is given by 


_ n - 27 Cn,~f 

2n ' 


( 1 . 22 ) 


(1.23) 


Hamiltonian quantities for fractional problems have been recently developed in 
the setting of one-dimensional solutions for fractional semilinear equations (—A)'''u'-|- 
f{w) = 0. The first reference we find a conserved Hamiltonian quantity for this type of 
non-local equations is the paper by Cabre and Sola-Morales [7] for the particular case 
7 = 1/2. The general case 7 € (0,1) was carried out by Cabre and Sire in [ 8 ]. On the 
one hand, in these two papers [ 8 , 7], the authors impose a nonlinearity coming from a 
double-well potential and look for layers (i.e., solutions that are monotone and have 
prescribed limits at infinity), and they are able to write a Hamiltonian quantity that 
is preserved. In addition, if one considers the same problem but on hyperbolic space, 
one finds that the geometry at infinity plays a role and the analogous Hamiltonian is 
only monotone (see [33]). 

On the other hand, if, instead, one looks for radial solutions for semilinear equa¬ 
tions, then Cabre and Sire in [ 8 ] and Frank, Lenzman and Silvestre in [31] have 
developed a monotonicity formula for the associated Hamiltonian. In the setting of 
radial solutions with an isolated singularity for the fractional Yamabe problem, our 
Theorem 1.2 states that, if one uses the metric (1.19) to rewrite the problem, then the 
associated Hamiltonian (1.21) is constant along trajectories. 


If one insists on performing an ODE-type study for the PDE problem (1.20), a 
possibility is to plot a phase portrait of the boundary values (at p = 0 ), while keeping 
in mind that the equation is defined on the whole extension. From this point of view, 
one can prove the existence of two critical points: the constant solutions rip = 0 and 
= I. Moreover, there exists an explicit homoclinic solution r^oo, whose precise 
expression will be given in (6.1); it corresponds to the n-sphere (1.5). 

The next step is to linearize the equation. As we will observe in Section 7, the 
classical Hardy inequality, rewritten in terms of the background metric go , decides the 
stability of the explicit solutions uq, vi and Voo- Stability issues for semilinear frac¬ 
tional Laplacian equations have received a lot of attention recently. Some references 
are: [5] for the half-Laplacian, [59] for extremal solutions with exponential nonlinear¬ 
ity, [26] for semilinear equations with Hardy potential. In the particular case of the 
fractional Lane-Emden equation, stability was considered in [22, 29], for 7 S (0,1) and 
[28] for 7 G (1,2). We believe that our methods, although still at their initial stage. 
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would provide tools for a unified approach for all 7 £ (O, §)■ 

Finally, we consider the linearization of equation (1.12) around the equilibrium 

Vi = 1 : 

P^°v = on R X 

and look at the projection over each eigenspace {Ek), k = 0,1,.. 

( 1 - 24 ) 

Although we will not provide a complete calculation of the spectrum, we can say the 
following: 

Theorem 1.3. For the projection k = 0, equation (1.24) has periodic solutions v{t) 
with period where Xj is the unique positive solution of (7.3). In addition, 

lim Lj = Li, 

7-».l 

SO we recover the classical case (1.18). 

We conjecture that is the minimal period of the radial periodic solutions for 
the nonlinear problem ( 1 . 12 ). 

Remark 1.4. We also give some motivation to show that the projection on the k- 
eigenspace of (1.24) does not have periodic solutions if k = 1,2,_ 

Theorem 1.3 gives the existence of periodic radial solutions for the linear problem. 
In addition, the existence of a conserved Hamiltonian hints that the original non-linear 
problem has periodic solutions too. Moreover, the period behaves well under the 
limit 7 —^ 1 . Using the Hamiltonian and the methods in [ 8 ] one can prove that periodic 
solutions also behave well under this limit. We leave this proof for the forthcoming 
paper [23]. 

The construction of Delaunay solutions allows for many further studies. For in¬ 
stance, as a consequence of their construction one obtains the non-uniqueness of the 
solutions for the fractional Yamabe problem in the positive curvature case, since it 
gives different conformal metrics on S^(T) x that have constant fractional curva¬ 
ture. This is well known in the scalar curvature case (see the lecture notes [61] for an 
excellent review, or the paper [60]). In addition, this gives some motivation to define 
a total fractional scalar curvature functional, which maximizes the standard fractional 
Yamabe quotient ([35]) across conformal classes. We hope to return to this problem 
elsewhere. 

From another point of view, Delaunay solutions can be used in gluing problems. 
Classical references are, for instance, [50, 53] for the scalar curvature, and [51, 52] for 
the construction of constant mean curvature surfaces with Delaunay ends. 

There is an alternative notion of fractional curvature and fractional perimeter de¬ 
fined from the singular integral definition of the fractional Laplacian which gives a 




different quantity than our Q®. [ 11 ] introduces the notion of nonlocal mean curvature 
for the boundary of a set in R" (see also the review [64]), and it has also received a lot 
of attention recently. Finding Delaunay-type surfaces with constant nonlocal mean 
curvature has just been accomplished in [ 6 ]. For a related nonlocal equation, but 
different than nonlocal mean curvature, the recent paper [ 20 ] establishes variationally 
the existence of Delaunay-type hypersurfaces. 

We finally comment that the negative fractional curvature case has not been ex¬ 
plored yet, except for the works [17, 16]. They consider singular solutions for the 
problem (—A)'^w -I- = 0 in a domain fl with zero Dirichlet condition on dfl. 

This setting is very different from the positive curvature case because the maximum 
principle is valid here. We also cite the work [56], where they consider singular so¬ 
lutions of AW = 0 in a domain with a nonlinear Neumann boundary condition 
d^W = f{x,W) -W on dn. 

The paper will be structured as follows: in Section 2 we will recall some standard 
background on the fractional Yamabe problem. In particular we present the equivalent 
formulation as an extension problem coming from scattering theory. In section 3 we 
will give a geometric interpretation of the problem. Next, in section 4 we will analyze 
the scattering equation to give a proof for theorem (I.l). That is, we will compute the 
Fourier symbol for the conformal fractional Laplacian. In section 5 we face the problem 
from an ODE-type point of view. This kind of study over the extension problem (1.20) 
gives us two equilibria and the existence of a Hamiltonian quantity conserved along 
the trajectories. Moreover we will find in Section 6 an explicit homoclinic solution, 
which corresponds to the n—sphere. Finally, in Section 7 we will perform a linear 
analysis close to the constant solutions which corresponds to the n-cylinder. 

2 Preliminaries 

The conformal Laplacian operator for a Riemannian metric 5 on a n-dimensional 
manifold M is defined as 

Lg = -Ag + CnRg, where c„ = (2.1) 

and Rg is the scalar curvature. The conformal Laplacian is a conformally covariant 

4 

operator, indeed, given and g two conformally related metrics with = w^-^g, 
then the operator Lg satishes 

LgAf) = w~^Lg{wf), 

for every / G C°°{M). In the case / = 1 we obtain the classical scalar curvature 
equation 

n+2 

-AgW + CnRgW = CnRg^W^-2 , 

which, in the flat case, is precisely equation (1.15). 
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The fractional Laplacian on R" is defined through Fourier transform as 


{—A)'^w = V 7 G M. 

Note that we use the Fourier transform defined by 

{(;(^) = (27r)“"/^ / w{x)e~^^'^ dx. 

JR" 

Let 7 G (0,1) and u G L°° fl in M", the fractional Laplacian in R" can also be 
defined by 




w{x + y) — w{x) 
|y|("+27) 


dy, 


where P.V. denotes the principal value, and the constant (see [48]) is given by 




Caffarelli-Silvestre introduced in [12] a different way to compute the fractional 
Laplacian in R" for 7 G (0,1). Take coordinates x G R”, y G R+. Let w be any 
smooth function defined on R" and consider the extension W : R" x R+ —>■ R solution 
of the following partial differential equation: 


div(y“VLF) =0, a: G R”, ?/ G 
W(x,0) = w(x), a: G R", 


( 2 . 2 ) 


where a = 1 — 2j. 
kernel 


Note that we can write W = *x w, where is the Poisson 


K'y{x,y) = c 


i\x\ 


(2.3) 


and c = c(n,"f) is a multiplicative constant which is chosen so that, for all j/ > 0 , 
JK, {x, y) dx = 1. In addition. 


(—A)^in = —dj lim y'^dyW, 

y^0+ 


for the constant 

^ 7r(-7) ■ 


(2.4) 


One can generalize this construction to the curved setting. Let be a smooth 

manifold of dimension n + 1 with smooth boundary dX = M". A defining function in 
X for the boundary M is a function p which satisfies: 


p > 0 in A, p = 0 on M and dp ^ 0 on M. (2.5) 

A Riemannian metric g~^ on X is conformally compact if (A,p) is a compact Rie- 
mannian manifold with boundary M for a defining function p and g = p'^g'^ ■ Any 
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conformally compact manifold {X,g^) carries a well-defined conformal structure [g\ 
on the boundary M, where g is the restriction of g\M- We call (M, [g]) the confor¬ 
mal infinity of the manifold X. We usually write these conformal changes on M as 

4 

gw = w g ioT a. positive smooth function w. 

A conformally compact manifold {X,g'^) is called conformally compact Einstein 
manifold if, in addition, the metric satisfies the Einstein equation RiCg+ = —ng'^, 
where Ric represents the Ricci tensor. One knows [37] that given a conformally com¬ 
pact Einstein manifold {X,g~^) and a representative g in [ 5 ] on the conformal infinity 
M”, there is an unique defining function p such that one can find coordinates on a 
tubular neighborhood M x (0, e) in X in which g'^ has the normal form 

g+= p-'^{dp^ + gp), (2.6) 

where gp is a family on M of metrics depending on the defining function and satisfying 
9p\m = 9- 

Let {X,g~^) be a conformally compact Einstein manifold with conformal infinity 
(M, [ 5 ]). It is well known from scattering theory [39, 49, 40] that, given w G C°°{M) 
and s £ C, if s(n — s) does not belong to the set of L^-eigenvalues of — Ag+ then the 
eigenvalue problem 

-Ag+U - sin-s)U = 0 in X, (2.7) 

has a unique solution of the form 

U = Wp^-^ + Wip\ W,WiGC^(X), W\p^o = w. (2.8) 

Taking a representative g of the conformal infinity (M, [g \) we can define a family of 
meromorphic pseudo-differential operators iS'(s), called scattering operators, as 

Sis)w = Wi\M. (2.9) 

The case that the order of the operator is an even integer was studied in [39]. More 
precisely, suppose that m £ N and m < ^ if n is even, and that (f )^ — nn? is not an 
eigenvalue of —Ag+, then S'(s) has a simple pole at s = ^ -I- m. Moreover, if 
denotes the conformally invariant GJMS-operator on M constructed in [38] then 

CmP^ = - Res,= ^ Sis), Cm = (-1)™(m - l)!]-\ 

where ReSs^sj, S'(s) denotes the residue at sq of the meromorphic family of operators 
Sis). In particular, if m = 1 we have the conformal Laplacian Pf = Lg from (2.1), 
and if m = 2, the Paneitz operator 

Pi ~ i~^9) P S{anRg + bnRicg)d H- - — Qg^ 

where Qg is the Q-curvature and a„, are dimensional constants ([57]). 

It is also possible to define conformally covariant fractional powers of the Laplacian 
in the case 7 ^ N. For the rest of the paper, we set 7 £ (O, ^) not an integer, s = f-I- 7 . 
In addition, assume that s(n — s) is not an L^-eigenvalue for —Ag+ and that the first 
eigenvalue Ai(— Ag+) > s(n — s). In this setting: 
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Definition 2.1. We define the conformally covariant fractional powers of the Lapla- 
cian as 

PS[g+,g] = d^S{^+j), where (2.10) 

As a pseudodifferential operator, its principal symbol coincides with the one of 
(—Ag)^. In the rest of the paper, once g~^ is fixed, we will use the simplified notation: 

Py ■= P~f[9^,9]- 

These operators satisfy the conformal property 

P9»/ = U;-^P9(^/), WfGC°°{M), (2.11) 

for a change of metric 

4 

gtn ■■= w > 0 . 

Definition 2 . 2 . We define the fractional order curvature as: 


Q% ■■= ^"(1)- 


Note that up to multiplicative constant, Qi is the classical scalar curvature and 
<52 is the so called Q-curvature. 

Remark 2.3. Using the previous definition we can express the conformal property 
( 2 . 11 ) as 

n + 2'y 

pg^w) =w^QS-. ( 2 . 12 ) 

Explicit formulas for P® are not known in general, however, Branson [4], gave an 
explicit formula for the conformal Laplacian on the standard sphere, i.e. 


pSs" 

7 


r(i? + 7 + i) 
r(p-7 + i)’ 


(2.13) 


where B = y — Agg„ + (^^)^- For example, 

pSs" _ _ A I ri(n-2) pQs^ _ /_ A , ^ n-1 ^ 

-“1 ~ ~r 4 ’ ^1/2 ~ \J V 2 / ■ 

From (2.13) we can compute the fractional curvature on the unit sphere as 

rf 2 i 4-A,) 

QT =Pr(^)= ( 2 - 14 ) 

It is proven in [14] (see also the more recent paper [13]) that the conformal frac¬ 
tional Laplacian is the Dirichlet-to-Neumann operator for an extension problem that 
generalizes ( 2 . 2 ): 
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Theorem 2.4. Lei 7 G (0,1) and {X,g~^) be a conformally compact Einstein manifold 
with conformal infinity {M,[g]). For any defining function p of M satisfying (2.6) in 
X, the scattering problem (2.7)-(2.8) is equivalent to 


- div(p“VVF) + E{p)W 
W 


= 0 m {X,g), 
= w on M, 


(2.15) 


where 

g = P^g+, W = p^-^U, s=t + 7, 0 = 1-27. 

and the derivatives in (2.15) are taken respect to the metric g. The lower order term 
is given by 

E{p) = -A,(pt)pt + ( 7 ^ - i) 
or written back in the metric g~^, 


E{p) = p-^-\-Xg+ - s{n - s))p^-\ (2.16) 


In addition, we have the following formula for the calculation of the conformal frac¬ 
tional Laplacian 


Pfiw = —dry lim p°'dpW, 


where dry is defined as 


22'>'-ir(7) 

7r(-7) 


Remark 2.5. If X is the hyperbolic space identified with the upper half space 

with the metric g^ = , then the conformal infinity is simply M = R" 

with the standard Euclidean metric Idccp and therefore, problem (2.15) is precisely 
the extension problem considered by Caffarelli-Silvestre (2.2). As a consequence, the 
conformal fractional Laplacian reduces to the standard fractional Laplacian without 
curvature terms, i.e., = (—A)^. 

Now we are going to choose a suitable defining function p*, in order to transform 
the problem (2.6) into one of pure divergence form. We follow [14, 13]: 


Theorem 2.6. Set^ G (0,1)- Let (X, g^) be a conformally compact Einstein manifold 
with conformal infinity [M, [gr]), and such that Ai(—Ag+) > tL —^ 2 ^ Assuming that p 
is a defining function satisfying (2.5), there exists another (positive) defining function 
p* on X, such that for the term E[p) defined in (2.16) we have 


E{p*) = 0. 


The asymptotic expansion of this new defining function is 


P* = P 1 + 


2 Q® 


(n - 27 ) 


P^" + 0(P^)). 
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In addition, the metric g* = {p*)'^g'^ satisfies g*\p^o = g and has asymptotic expansion 
g* = [dp*f[l + 0[[p*f^)] + g[l + 0{{p*n 


The scattering problem (2.8)-(2.7) is equivalent to the following one. 


-div((p*)“VVF) 

W 


= 0 in iX,g*), 
= w on M, 


where the derivatives are taken with respect to the metric g* and W = {p*y "C/. 
Moreover 

P^w = -d^ W + wQy (2.17) 

The fractional Yamabe problem is, for 7 £ (0)§)) to find a new metric g^ = 

4 

y]n- 2 ~, g OH M conformal to g, with constant fractional curvature Q®". Using the con¬ 
formal property (2.12) the Yamabe problem is equivalent to find w a strictly positive 
smooth function on M satisfying 

n-\-2y 

P®(w) = cw"- 2 t , w > 0 . (2-18) 

In this paper we are interested in the positive curvature case, and the constant c = 
will be normalized as in Proposition 2.7 below. 

Thanks to Theorem 2.4, (2.18) is equivalent to the existence of a strictly positive 
C°° solution for extension problem: 


- div(p“VlU) -k E{p)W = 0 in {X, g), 

W = w on M, 

71 + 27 

—dy lim p^dpW = Cn on M 

’ p^O ^ ' 


(2.19) 


Using the special defining function from Theorem 2.6, the fractional Yamabe problem 
(2.19) may be rewritten as 


r -div((p*)“VlU) = 0 in(Y,g*), 

J W = wonM, ( 2 . 20 ) 

—d.y lim^(p*)“clp. lU -I- wQ® = on M. 

Indeed we only need to rewrite the equation for the Yamabe problem (2.18) using the 
expression of P® from (2.17). Without danger of confusion, note that in general the 
solutions W (2.19) and (2.20) are different, but in the sequel they will be denoted by 
the same letter. 

Proposition 2.7. The fractional curvature of the cylindrical metric g^.,^ = wi \dx\'^ 
for the conformal change 

wi{x) = |a:| 2 -^, ( 2 . 21 ) 
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is the constant 


2 


Proof. The value is calculated using the conformal property (2.12), as follows: 

Q®”" = (wi) = (-A)'^(r(;i) =: 

The last equality follows from the calculation of the fractional Laplacian of a power 
function; it can be found in [34, 58]. □ 


3 Geometric setting 


We give now the natural geometric interpretation of problem (1.1) and the extension 
formulation (1.3). Thanks to Remark 2.5 and Theorem 2.4, the initial extension 
problem (1.3) is the scattering equation (2.7) in hyperbolic space, denoted by Xi = 
]HI”+^, with the metric ■ Our point if view is to use the metric go 

from (1.10) as the representative of the conformal infinity, thus we need to rewrite the 
hyperbolic metric in a different normal form 


5 "^ = 


dp^ + gp 


with gp|p=o = go, 


(3.1) 


for a suitable defining function p. Let us introduce some notation now. The conformal 
infinity (with an isolated singularity) is Mi = IR"\{0}, which in polar coordinates can 
be represented as Mi = 1R+ x and \dx\'^ = or using this change of 

variable r = e“*, the Euclidean metric may be written as 


\dx\'^ = e ^‘[dt^+gs"-i] =: e ^‘go- 


(3.2) 


We consider now several models for hyperbolic space, identified with the Riemannian 
version of AdS space-time. These models are well known in cosmology since they 
provide the simplest background for the study of thermodynamically stable black holes 
(see [65, 45], for instance, or the survey paper [15]). Thus we write the hyperbolic 
metric as 

g^ = dcr^ -I- cosh^ a dt^ -\- sinh^ crggn-i, (3.3) 

where t € K., cr G (0, oo) 6 G Using the change of variable R = sinh cr, 

9^ = ^ ^2 -|- (1 -k d?^) dt^ + R^gsn-i. 

This metric can be written in the normal form (3.1) as 


= P 


<iy + (i + <i)b(> + (i-4)’ 


(3.4) 
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for p € (0, 2), t G M, 0 £ §" Here we have used the relations 

p = 2e"‘" and 1 + (3.5) 

Let g = be a compactification of g^. Note that the apparent singularity at p = 2 
in the metric (3.4) is just a consequence of the polar coordinate parametrization and 
thus the metric is smooth across this point. 

We define now X 2 = (0,2) x §^(L) x with coordinates p G (0,2), t G 

S^(L), 9 G and the same metric given by (3.4). The conformal infinity {p = 0} 

is M 2 = S^(L) X S"“^, with the metric given by po = dt^ + 

Note that (Xi,pj|(„+i) is a covering of (X 2 ,p+). Indeed, X 2 is the quotient X 2 = 
IHI"+^/Z Ri K." X S^(L) with Z the group generated by the translations, if we make the 
t variable periodic. As a consequence, also (Mi, \ dx\^) is a covering of {M 2 , go) after 
the conformal change (3.2). 

Summarizing, we denote X = (0, 2) x R x S"“^ and M = R x and recall that 
the metric g = p^g'^ is given by 

p = dp2 + ^l + ^) pgn-i, and po = 5 |m = +5 §"-i- (3.6) 

Equality (3.2) shows that the metrics jfixp and po are conformally related and therefore 
using (1.6), we can write any conformal change of metric on M as 

g^ ■= w^^\dx\'^ = go- (3.7) 

Our aim is to to find radial (in the variable |x|), positive solutions for (1.3) with an 
isolated singularity at the origin. Using po as background metric on M, and writing 
the conformal change of metric in terms of v as (3.7), this is equivalent to look for 
positive solutions v = v{t) to ( 1 . 12 ) with 0 < ci < u < C 2 , and we hope to find those 
that are periodic in t. 

Finally, we check that the background metric po given in (3.6) has constant frac¬ 
tional curvature (5®° = Cn,-y This is true because of the definition of given in 
Proposition 2.7, and the conformal equivalence given in (3.2). Thus, by construction, 
the trivial change = 1 is a solution to ( 1 . 12 ). 


4 The conformal fractional Laplacian on M x S” 

In this section we present the proof of Theorem 1.1, i.e, the calculation of the Fourier 
symbol for the conformal fractional Laplacian on R x This computation is based 

on the analysis of the scattering equation given in (2.7)-(2.8) for the extension metric 
(3.4). We recall that the scattering operator is defined as 

psy: = S{s)w = Wi\p=o, (4.1) 
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and s = § + 7- The main step in the proof is to reduce (2.7) to an ODE that can 
be explicitly solved. Note that this idea of studying the scattering problem on certain 
Lorentzian models has been long used in Physics papers, but in general it is very hard 
to obtain explicit expressions for the solution and the majority of the existing results 
are numeric (see, for example, [19]). 

For the calculations below it is better to use the hyperbolic metric given in the 
coordinates (3.3). Then the conformal infinity corresponds to the value {a = +oo}. 
The scattering equation (2.7) can be written in terms of the variables a G (0, oo), 
t gR and 9 G as 

daaU + Q{(7)daU + cosh~‘^{a)dttU + sinh“^(cr)Asn-il7 + (x “ ^ = 0; (4-2) 

where U = U{a,t,9), and 


Q(o-) = 


With the change of variable 


dcr (cosh cr sinh" ^ a) 
cosha sinh"“^ a 

z = tanh(cr). 


equation (4.2) reads: 

(1 - +{!i^-z){l- z^)d,U + (1 - z^)dttU 

+ (^-1) As„-it/+(2^-7')t/ = 0. 


(4.3) 


(4.4) 


We compute the projection of equation (4.4) over each eigenspace of Agn-i. Given 
fc G N, let Uk{z,t) be the projection of U over the eigenspace (Ek) associated to the 
eigenvalue = —k(k + n — 2). Each Uk satisfies the following equation: 


(1 ~ z^)dzzUk + ~ -s) dzUk + dttUk + fJ-k-^Uk H— fzr^Uk — 0. 


Taking the Fourier transform (1.13) in the variable t we obtain 


{l-Z^)dzzUk+{^-z)dzUk+ Pik-^ 




-e 


Uk = 0. 


(4.5) 


(4.6) 


Fixed k and we know that 

Uk = wi{^)(fl{z), 

where Lp := V^|(z) is the solution of the following ODE problem: 

f (1 — z^)dzz'P + ~ z) dz'p + + \_Ji -7^ = 0, 

I has the expansion (2.8) with w = 1 near the conformal infinity z = 1, 
[ is regular at z = 0. 


(4.7) 


(4.8) 
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This ODE has only regular singular points z. The first equation in (4.8) can be 
explicitly solved, 


ip{z) =A{1 - Hyperg(a,6;c; z^) 

+B(1 - Hyperg(a,6; 

for any real constants A,B, where 


(4.9) 


a = ^ 

“ 2 ' 2 


1 , t /(2 , •? 


+ * 2 ’ 


2 ~'~ 2 ~ 2 


C = 1 + v'd -l) 2 -^fc, 


• d = ^ , 2 


b = ^ , 

" 2 ' 2 


1 _ t/(2 j_ -i 

2 2 

1 


+ * 2 ’ 


— If 


• c = 1 - Vd -1)2 - Mfc> 

and Hyperg denotes the standard hypergeometric function introduced in Lemma 4.3. 
Note that we can write ^ instead of |^| in the arguments of the hypergeometric functions 
because a = b, a = b and property (4.19) given in Lemma 4.3. 

The regularity at the origin z = 0 implies i? = 0 in (4.9). Moreover, property 
(4.18) from Lemma 4.3 makes it possible to rewrite (p as 

ip{z) = A a{l + z)'^~^ (1 — z)'^~^ z^ 2+\/(~2 1) 

• Hyperg(a, 6 ; a + 6 — c + 1; 1 — z^) 

+ /3(1 + z)?+^ (1 - 

• Hyperg(c — a,c — 6 ;c — a — 6 +l;l — z^)] , 
r(i+/(f^^Ipd^)r(7) 


(4.10) 


where 




^(2 1)^ 


— it; I r I i + ? + 


•\/( 2 . g 


(4.11) 


/3 = 






— 7 1 

^+2 + 


.5' 


The constant coefficient A will be fixed from the second statement in (4.8). From the 
definition of the scattering operator in (4.1), p must have the asymptotic expansion 
near p = 0 

p{p) = p"-7l + ...) + p%SHs)l + ...), (4.12) 
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where S^{s) is the projection of the scattering operator S{s) over the eigenspace (Ek). 
We now use the changes of variable (4.3) and (3.5), obtaining 

z = tanh(cr) = = 1 - . (4.13) 

Therefore, substituting (4.13) into (4.10) we can express as a function on p as follows 

ip{p) ~ A [ap'2 Hyperg(a, b; a + b — c + 1; p^) 

+ Pp^^^ Hyperg(c — a, c — 6; c — a — 6+1; p^)] , as p —0. 

Using property (4.17) from Lemma 4.3 below, we have that near the conformal infinity, 

ip{p) ~ A [ap'^~^ + Pp^'^^ + ...]. (4-14) 

Comparing (4.14) with the expansion of ip given in (4.12), we have 

A = a-\ (4.15) 

and _ 

S^is)=Pa-\ (4.16) 

Recalling the definition of the conformal fractional Laplacian given in (2.10), and tak¬ 
ing into account (4.7), we can assert that the Fourier symbol 0(((O for the projection 
P!^ of the conformal fractional Laplacian satisfies 

®7(^) = I^22'^^(s). 

From here we can calculate the value of this symbol and obtain (1.14); just take (4.16) 
into account and property (4.20) from Lemma 4.4. This completes the proof of The¬ 
orem 1.1. 


Remark 4.1. When 7 = to, an integer, we recover the principal symbol for the GJMS 
operators P^ ■ Indeed, from Theorem 1.1 we have that for any dimension n > 2m, 
the Fourier symbol of P^ is given by 


0^(0 = 2 


2m 


ir(i+f + dlEEEE+to 


|r(i - f + 


\/*^2 , « ■ 




= 2 


2m 





= 4-(TO,n,fe,e,e^-,r™-') + e 


t-2m 


where we have used the property (4.21) of the Gamma function given in Lemma 4.4. 
Note that 4' is a polynomial function on ^ of degree less than 2m. 
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For instance, for the classical case m = 1, 

so we recover the usual conformal Laplacian Pf° given by 


P^{v) = -v + 


(n-2f 


fJ-k 


V, k = 0,l,..., 


that is the operator appearing in (1.16) when applied to radial functions. 

This proof also allows us to explicitly calculate the special defining function p* 
from Theorem 2.6: 

Corolary 4.2. We have 


{p*r-^ = a-^ ( 4 ^ 


"■"Hyperg(f-2,|-l;§,(|^) ), 


where a is the constant from (4.11). As a consequence, p* € (0, pg) where we have 
defined (Po)"~* = a~^. 

Proof. From the proof of Theorem 2.6, which corresponds to Lemma 4.5 in [14], we 
know that 

p* = ((^O)^i^(z), 

where <p is the solution of (4.8). Thus from formula (4.9) for i? = 0 and the relation 
between 2 and p from (4.13) we arrive at the desired conclusion. The behavior when 
p —2 can be calculated directly from (4.9) and, as a consequence, (po)"~* = 'F(O) = 
a-i. □ 

Lemma 4.3. [2, 63] Let z € C. The hypergeometric function is defined for jzj < 1 

by the power series 


Hyperg(a, 6 ;c; z) = '^ 


{a)n{b)n z" _ _ 

(c)„ n! F(a)F( 6 ) 


r(c) r(a + n)T{b + n) z'' 


F(c + n) n! ' 


It is undefined (or infinite) if c equals a non-positive integer. Some properties are 

1. The hypergeometric function evaluated at z = 0 satisfies 

Hyperg(a + j, 6 - j; c; 0) = 1; j = ±1, ±2,... (4.17) 

2. If |ar 5 (l — z)\ < tt, then 

rfcirfc CL 

Hyperg(a, b; c; z) = ) ^ ^ Hyperg (a, 6; a + 6 - c + 1; 1 - z) 

i (c — ajl (c — b) 

+(1 ~ —~ Hyperg(c — a, c — 6; c — a — 6 + 1; 1 — z). 


r(a)r(6) 


(4.18) 
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3 . 


The hypergeometric function is symmetric with respect to first and second argu¬ 
ments, i.e 

Hyperg(a, 6 ; c; z) = Hyperg(&, a; c; z). (4.19) 


Lemma 4.4. [2, 63] Let z G C. Some properties of the Gamma function r(z) are 


r(z) = r(z), 

r(z + 1) = zT{z), 

T{z)T {z+\)= 2i-2^v^r(2z). 

Let tjj{z) denote the Digamma function defined by 

,, , dlnViz) r(z) 

= r(iT' 

This function has the expansion 

OO 

i,{z) = V'(l) + ^ . 

m=0 


Let B{zi,Z 2 ) denote the Beta function defined by 


B{zi,Z2) 


r(^i)r(z2) ^ 

r(zi + Z2)' 


(4.20) 

(4.21) 

(4.22) 


(4.23) 


If Z 2 is a fixed number and zi > Q is big enough, then this function behaves 


B{zi,Z2) ^T{z2){zi) 


(4.24) 


We end this section with a remark on the classical fractional Hardy inequality. On 
Euclidean space (K”, \dx\'^), it is well known that, for all w G C§°(R^) and 7 G (O, §), 


CH 


,|27 


dx < 






f \{-A)2w\'^dx 
Jr-^ 


w{—Ayw dx. 


(4.25) 


Moreover, the constant ch is sharp (although it is not achieved) and its value is given 

by 

Ch — Cn,7; 

which is the constant in Proposition 2.7. This is not a coincidence, since the functions 
that are used in the proof of the sharpness statement are suitable approximations of 
(2.21). This constant was first calculated in [42], but there have been many references 
[ 66 , 3, 32], for instance. 

A natural geometric context for the fractional Hardy inequality is obtained by 
taking go as a background metric, and using the changes (1.6) and (1.10). Indeed, 
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using the conformal relation given by expression (1.11), we conclude that (4.25) is 
equivalent to the following: 

c „,7 / dvolgfj < / v{P^°v)dvolgg, (4.26) 

for every v G C((°(IR x S"“^). 


5 ODE-type analysis 


In this section we fix 7 G (0,1). As we have explained, the fractional Yamabe problem 
with an isolated singularity at the origin is equivalent to the extension problem (1.3). 
We look for radial solutions of the form (1.6). Based on our previous study, it is 
equivalent to consider solutions v = v{t) of the extension problem ( 1 . 20 ), for the 
metric (1.19). In this section we perform an ODE-type analysis for the PDE problem 
( 1 . 20 ). 

Firstly we calculate 


1 


diVgip^y-gV) =Y^ 


\/M 

= (p“e(p)5pE) -I- 


(5.1) 


P , dttV + P Agn-iE, 

( 1 + 4)2 ( 1 - 4)2 


where 

«(»>) = (1 + ^) (1 - i) 

Using the expression given in (2.16), 


E{p) 


n—1+a 
4 ^ 


n — 2 -I- n4 

MIF?)' 


(5.2) 


Remark 5.1. Let V be the (unique) solution of (1.20). If v does not depend on the 
spherical variable 6 G then V does not either. Analogously, if v is independent 

on t and 9, then V is just a function of p. The proof is a straightforward computation 
using that the variables in (5.1) are separated. 

As a consequence of the previous remark, it is natural to look for solutions V of 
(1.20) that only depend on p and t, i.e. solutions of 


-^5p(p“(e(p)5pU) 


(1 + f)^ 


dttV + Eip)V = 0 for pG (0,2),t G 
V = V on {p = 0}, 


—d-y lim p°'dpV = Cn ■yV on {p = 0}. 

p —>-0 ’ 


(5.3) 
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Now we take the special defining function p* given in Theorem 2.6, whose explicit 
expression is given in Corollary 4.2. Then we can rewrite the original problem (1.20) 
in g*, defined on the extension X* = {p G (0, Pq), t G K, 0 G as 

-div,.((p*)“Vg.y) 

c 

—d^ hm {p*)°'dp»V + Cn^-yV 


0 in(X*,p*), 

V on {p* = 0}, 

n + 2~Y 

c„, 7 f on {p* = 0 }, 


(5.4) 


where g* = for p* = p*(p). 

Note that Proposition 2.7 calculates the value = Cn,-y The advantage of (5.4) 
over the original ( 1 . 20 ) is that it is a pure divergence elliptic problem and has nicer 
analytical properties. 

Next, if we look for radial solutions (that depend only on t and p*), then the 
extension problem (5.4) reduces to: 


^ dp, {{p*Te*{p)dp,V) + ^ duV = 0 for t G ]R,p* G (0,pS), 


s*(p) 




= fo on {p* = 0 }, 


(5.5) 


-dj{p*)°'dp,V + Cn,^yV =Cn,jV”-'^-> On {p* = 0}, 


where 

e*{p)={^y e{p). 

Summarizing, we will concentrate in problems (5.3) and (5.5). In some sense (5.5) 
is closer to the local equation (1.16) and shares many of its properties. For instance, it 
has two critical points: uq = 0 and ui = 1 , since these are the only constant solutions 

n+2~f 

of the boundary condition v = on p* = 0. Moreover, by uniqueness of the 

solution and Remark 5.1, the only critical points in the extension are simply Vq = 0 
and hi = 1 . 

Remark 5.2. The calculation of the critical points uq = 0 and = 1 also holds for 
any 7 G (0,since the corresponding extension problem shares many similarities 
with (5.5) (c.f. [14, 13, 67]). 

The linearization at Vi = 1 will be considered in Section 7. 


5.1 A conserved Hamiltonian 


Here we give the proof of Theorem 1.2. The idea comes from [ 8 ], where they consider 
layer solutions for semilinear equations with fractional Laplacian and a double-well 
potential. Multiply the first equation in (5.3) by e{p)dtV, and integrate with respect 
to p G (0, 2), obtaining 


- dp{p^e{p)dpV)dtVdp- p^ei{p)duVdtVdp+ p^e 2 {p)VdtV dp = Q, 
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where we have defined e, ei, 62 as in (1.22). We realize that dttVdtV = ^dt{{dtV)^) 
and VdtV = thus integrating by parts in the first term above we get 

p°‘e{p)dpVdtpV dp + {p°-e{p)dpVdtV) |p=o 

J dp^ J p°‘e2{p)V'^ df^ = 0 . 

Here we have used the regularity of H at p = 2. Again we note that dpVdtpV = 
^dt{{dpV)^) and using the boundary condition, i.e., the third equation in (5.3), we 
have 



^dt p‘^e{p){dpVf df^ - \dt p‘^ex{p){dtVf dp 

* (i 


(5.6) 


P°‘^2{p)V^dp] = 

d'y 


Define 

G{v) = , 

where the constant is defined in (1.23). In this way, we have from (5.6) that 

^dt [ {p‘'e{p){dpVf - p°‘ei{p){dtVf + p“e 2 (p)H^} dp - dt{G{v)) = 0. 

Jo 

So we can conclude that the Hamiltonian 

-=1 ( {e{p){dpVf - ei{p){dtVf + e 2 (p)H^} dp - G(v), 

Jo 

is constant respect to t. This concludes the proof of Theorem 1.2. 

Remark 5.3. One can rewrite the Hamiltonian in terms of the defining function p*. 
For this, we may follow similar computations as above but starting with equation 
(5.5). Indeed, let H be a solution of (5.4), then the new Hamiltonian quantity 


Hyt) : 


fPo 

+ 5 / iP*r {eUp){dtVr - e*{p){dp.vr} dp* 

Jo 


(5.7) 


is constant respect to t. Here 

e*{p) = [^)\ elip) = ay e,{p). 


This quantity H* is the natural generalization of (1.17). 
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Now we observe that in the local case, the Hamiltonian (1.17) is a convex function 
in the domain we are interested, thus its level sets are well defined closed trajectories 
around the equilibrium ui = 1. We would like to have the analogous result for the 
Hamiltonian quantity H* from (5.7). This is a very interesting open question that we 
conjecture to be true. In any case, the second variation for H* near this equilibrium 
is: 


de^ 


H^iVi+eV) 

e=0 


n -27 ^ 


+ i 


fPo 


(P*)“^ {eiipKdtVf - e{p){dp.V)^} dp*. 


6 The homoclinic solution 

For this section we will take 7 £ (O, ^), since it does not depend on the extension 
problem (1.3). It is clear that the standard bubble (1.5) is a solution of equation (l.I) 
that has a removable singularity at the origin. Note that, because of our choice of 
the constant we need to normalize it by a positive multiplicative constant. We 
prove here that, on the boundary phase portrait, the equilibrium = I stays always 
bounded by this homoclinic solution in a boundary phase portrait. More precisely: 

Proposition 6.1. The positive function 


n.-2-, ( r(?-7)\“^^ 

?^oo(t) = C'(cosht) 2 ^ with C = + 7 ) j >l = r'i, (6.1) 

is a smooth solution of the fractional Yamabe problem (1.20). The value of Cn,~j is 
given in Proposition 2.1. 

Proof. The canonical metric on the sphere, rescaled by a constant, maybe written as 
gc = = [C{cosht)~^^]'^ go. 

We choose C such that the fractional curvature of the standard sphere is normalized 
to 

Qf = Cn,j. (6.2) 

Now we use the conformal property (2.12) for the operator P®®”: 


P9^"(C) = . (6.3) 

One checks that the fractional curvature is homogeneous of order 7 under rescaling of 
the metric. Indeed, because of (6.3) and the linearity of the operator P^ 

QSC = c- — p3^^ (C) = " (I) = C" —. (6.4) 
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Comparing equalities (6.2) and (6.4), together with the value of the curvature on the 
standard sphere (2.14) we find the precise value of C as claimed in (6.1). 

Next, let us check that the value of the constant C is larger than one. Because of 
Proposition 2.7 we have to test that 

o2.(^ r(i(§ + 7)) \^ r(f-7) 

lr(i(i-7)); r(f + 7) • 


Using the property (4.22) of the Gamma function, given in Lemma 4.4, we only need 
to verify that 


X(n,7) 


r(Mi + 7))r(^(i-7) + ^) 

r(Ki-7))r(i(i + 7) + ^) 


Thanks to Lemma 6.2 below, it is enough to see that 


X(n, 0) < 1 Vn, 


which holds trivially. 

Lemma 6.2. The function X(n,"f) defined as follows 


□ 


X(n, 7 ) := 


r(Ki + 7))r(i(i-7) + ^) 


r(i(i-7))r(i(§ + 7) + i)’ 

is increasing in n, and decreasing in 7 . 

Proof. If we denote ip^z) the Digamma function from Lemma 4.4, we can use the 
expansion (4.23) to study the growth of the function X{n,j) with respect to n and 7 . 
First, 

|-(logX(n, 7 )) = i (V^(f + 2 ) + V;(f - 7 + i) - _ 2 ) _ V ;(2 + 2 + 1 )) 


= iE 


m+5 + i 


m—0 


[(m+f)=-^J[(m+f + i)^-^J 


> 0. 


and 

A 

dj 


(iogx(n,7)) = i im + + m -^)-m + i + 1)) 


= -^E 


m—O 


(m+2 + l)(n^+2) + ^ 


(6.5) 


< 0 . 


□ 
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7 Linear analysis 


Let us say a few words about stability. Let z;* be a solution of (1-12). The correspond¬ 
ing linearized equation is 


P^°v = c 


i-i 

, !i±27„,»-2T 


We say that u, is a stable solution of (1-12) if 


dvolgg — c. 


! M 


n +27 




4-j, 

■= 2 y, 


M 


' dvolgg > 0 


for all veC^{M). (7.1) 


We observe here that the equilibrium ui = 1 is not a stable solution for (1.12) just 
by comparing the constant appearing in (7.1) and in the Hardy inequality (4.26). In 
addition, one easily checks that the equilibrium solution uq = 0 is stable. 

But it is more interesting to look at the explicit solution Voo given in (6.1). It 
follows from the Hardy inequality (4.26) that this explicit solution is not stable. The 
kernel of the linearization at Uoo is calculated in [ 21 ], where they show that, although 
non-trivial, is non-degenerate, i.e., is generated by translations and dilations of the 
standard bubble. 


Let us look more closely at the spectrum of the operator It is well known 

that is self-adjoint ([39]), and then we can compute its first eigenvalue through 
the Rayleigh quotient. Thus we minimize 

lM^P7°'>^dvolgg 
vaC^iM) Jj^v'^dvolga 

where M = R x We can apply Theorem 4.2 and Corollary 4.3 in [35] (or the 

Hardy inequality (4.26)) to conclude that P®“ is positive-definite. Moreover, the first 
eigenspace is of dimension one. 


Now we consider the linear analysis around the equilibrium solution z;i = 1. In 
order to motivate our results, let us explain what happens in the local case 7=1 for 
the linearization (see [50, 54, 47]). In these papers the authors actually characterize 
the spectrum for the linearization of the equation 


Pf°v = 


- 

4 I 


given by (after projection over each eigenspace (Pfc), fc = 0,1,...) 


—V — [n — 2 -I- ^k]v = 0. 

Note that this equation has periodic solutions only for fc = 0, of period Li = for 
A° = n — 2. Thus we recover (1.18). For the rest of fc = 1,..., the corresponding 
A^' = n — 2 -I- < 0, so we do not get periodic solutions. 
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The linearization of equation (1.12) around the equilibrium tii = 1 is given by 


pgoy — f, 


(7.2) 


Here we will calculate the period of solutions for this linearized problem (for the 
projection fc = 0), as stated in Theorem 1.3, by the method of separation of variables. 
We also conjecture that there are not periodic solutions for the linearized problem 
(1.24) for the rest of fc = 1,..., as it happens in the classical clase. 

Therefore, we consider the projection of equation (4.2) over each eigenspace (i?fc), 
fc = 0,1,.... Let 

Uk{z,t)=T{t)Z{z), 

be a solution of (4.5). Then 






Z'jz) 

Z(z) 



T"(t) ,, 

Tit) 


for a constant := A^( 7 ) £ K. We are only interested in the case A > 0, which is the 
one that leads to periodic solutions in the variable t. The period would be calculated 
from ;= 

Note that the equation for Z(z) is simply (4.6) with replaced by A^. From the 
discussion in Section 4, in particular (4.10), (4.15) and (4.16) we have that 


Ziz) =(i + 2 )t-3(i - 

Hyperg(a, b; a + b — c + 1;1 — z'^) 

+tt;(l + z)T+i(l - 

Hyperg(c — a,c — 6;c—a — 5+l;l — z^), 


where 


a = 

b = 
c = 


K = 


-Mfc 


r(-7) 


1 ^ 2 ^ ^ 2 

L , = 

■ vw 

1 ' 2 

i 2 I 



.,1,7, 

^2 ' 2 ' 2 

' 2 ■ 

1 7 \/( 



r(7) 


We use the change of variable (4.13) to analyze the asymptotic behavior of Z near the 
conformal infinity p = 0 

z - pt-T" + KP?+T'. 
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From the definition of the scattering operator (2.8), (2.9), and the definition of the 
conformal fractional Laplacian we have that 




r(i + i + + 41 .) 




Imposing the boundary condition (7.2) and the value of Cn,-^ given in (2.7), the un¬ 
known A^’ must be a solution of 


r(4 


xAFi)^ 

2 


r(^-i + 


y/ 

2 


0 


2 "T 2 


n + 2-f\T (§ + 7 ))! 

”-2T'|r(i(i-7))| 


(7.3) 


Note that for the canonical projection fc = 0, equality (7.3) simplifies to 


r(| + 2 + 

r(f-i + ^*) 


n + 2 -i\T (i + 7 ))! 


(7.4) 


This equation (7.4) lets us recover the value of A° for the classical case 7 = 1 . Indeed, 
using property (4.21) we get A° = n — 2 and we recover (1.18), where Lx := T°(l). 

Going back to equation (7.3) we can assert that the value of A^ can not be zero 
and it is unique for each k. Indeed if A = 0 we get a contradiction, and if A > 0 we 
may proceed as follows. Define 


Fm 


|^(afc-^,a^)p 


|r(Kf-7))r 


where 


i + 


2 


1 2 




and (3 = 


y/ 


Note that equation (7.3) is written as F(/3) = 1, for some /3 > 0. We derive this 
expression with respect to (3, 


[log Fm' 


- . .',2 3[V'(afc + /3i) - ^{ak + Pi)]- 

n + 2 'y 1^(212 

|r(Kf-7))r 
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Here S represents the imaginary part of a complex number and the Digamma 
function from Lemma 4.4. We can use the expansion (4.23) to arrive at 


(logF(/3))' = c^ 

m—0 


7/3^2m+l + y(T 

1 

7 

-+2+^ ' 2 



for some positive constant c. Therefore F(P) is an increasing function of /3. 
Next, note that 


lim F(/3) = + 00 , 

^—>• + 00 


for all A: = 0,1,.... This follows easily writing 




and the asymptotic behavior for the Beta function (4.24) from Lemma 4.4. 
Now we look at the projection k = 0. One immediately calculates 


F(0) 


n — 27 
n + 27 


< 1, 


so there exists (and it is unique) a solution A° = A°( 7 ) > 0 for the equation F{(3) = 1. 
From the proof one also gets that 


lim A°( 7 ) =n-2. 


This concludes the proof of Theorem 1.3. 


We believe that, as in the classical case F{f3) = 1 does not have any positive 
solution for k = 1,2,.... This is a well supported conjecture that only depends on 
making more rigorous some numerical analysis. In order to motivate this conjecture, 
let us try to show that /^ > 1 for /c = 1 , 2 ,..., where we have defined 


i^(0) 


{n - 2j)\T{ak)\'^ 

|r(Hi-7))| 

|2 

(n + 27 ) r(afc)p 

r(Hi + 7))| 

|2 


= : fk- 


Using the same ideas as above, one checks that fk is an increasing function of fc, and 
it is enough to show that 


(n-27)|r(i + i + f)r(f-i) 
(, + 27)|r(i-i + f)r(f + i) 


n — 2"f 
n + 27 


X(n, 7 ) ^ 


> 1, 


where X{n, 7 ) is defined in Lemma 6.2. We have numerically observed that fi = fiij) 
is an increasing function in 7 . Since for 7 = 0 we already have that /i(0) = 1, we 
would conclude that /fc > /i > 1 , as desired. 
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